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Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of 
two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalci- 
fications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7-0.8 mm. Current VA 
systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane 
are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is 
about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry 
degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D 
point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric 
mean filters. To assess the filters' performance on the restored images, we use an image quality index that accounts for correlation 
loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener 
filter is 0.9 (when the index is 1 this indicates perfect restoration). This filter yielded the best result in comparison with the other 
ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 
0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted 
in the diff'erence-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient infor- 
mation to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information 
to spot the wires in the images. It is demonstrated that the phase of the theoretical and the experimental PSFs are dissimilar. This 
fact prevents VA image restoration with the current theoretical PSF. This study is a preliminary step towards understanding the 
restoration of VA images through the application of deconvolution filters. 
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1. Introduction 

Recent practices of medicine strongly rely on the visual- 
ization of inner organs in order to study the relationship of 
anatomic structure and physiological process of living beings. 
In this context, vibro-acoustography (VA) |1| has become a 
medical imaging tool for investigating calcifications in human 
arteries EllSIll, microcalcifications in the breast (IIIH, liver 
tumors |7|, bone structure |8|, among others applications 

The VA technique employs two co-focused ultrasound beams 
in megahertz range with slightly different frequencies, typically 
in the kilohertz range, to scan biological tissues (see Fig. [T] for 
the system description). The VA image formation process is 
related to the nonlinear interaction of the primary ultrasound 
waves (incident and scattered), which gives rise to a secondary 
acoustic wave at the difference-frequency L10,,11J. The gen- 
erated diff'erence-frequency wave is detected by a hydrophone 
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(typical bandwidth 1-70 kHz). The acquired signal is processed 
and used to form VA images. 

The interaction of the primary incident ultrasound waves 
gives rise to a highly collimated difference-frequency beam, 
known as the parametric array effect 1 12 |. It has been recently 
shown that the difference-frequency signal detected in the for- 
ward scattering direction is mostly due to the nonlinear inter- 
action of the primary scattered waves |11|. Hence, the para- 
metric array may behave as a background signal present in the 
diff'erence-frequency wave detected to form VA images. 

The VA method can be extended to generate several images 
at once using a multifrequency arrangement also known as mul- 
tifrequency VA |[T3l . In this configuration, N ultrasound beams 
(more than two) are produced and they overlap in the system 
focal region. Due to the nonlinear interaction of the ultrasound 
waves, N{N - l)/2 difference-frequency components arise in 
the medium and they are used to form multiple images with 
just one scanning process. 

The spatial resolution of VA described in terms of the sys- 
tem point spread function (PSF) is determined by the trans- 
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Figure 1: Description of a typical VA imaging system. The central- and the 
difference-frequency are ojq and Aoj, respectively. 

ducer geometrical properties and the ultrasound driving fre- 
quencies. Beamforming has a great impact on the system res- 
olution. Some schemes to generate two or more ultrasound 
beams for VA have been proposed based on a confocal trans- 
ducer lfT4ll . x-focal (crossing beams) ifTsTl , sector array (161, 
linear array |[T6j[T3, and reconfigurable array fTE\. The VA 
resolution does not to correspond to the difference-frequency 
beam waist due to the parametric array effect. Indeed, the sys- 
tem PSF has a narrower lateral width than that of a parametric 
beam produced by the VA confocal transducer. To verify this, 
we have computed the parametric beam with the theory pre- 
sented in Ref. 1 19] and found that this beam does not fit the 
PSF waist shown in Ref. flAj . However, for the sake of brevity, 
these results will not be included in this article. 

The study presented here is limited to VA systems based 
on the confocal arrangement. Specifically, a confocal VA 
system composed by a two-element confocal transducer with 
44 mm-diameter, focused at 70 mm, and operating with fre- 
quencies around 3.1 MHz, achieves a spatial resolution (-6dB 
beam width) of approximately 0.8 x 0.8 nmi in the focal plane 
and 16 mm axially. With this resolution only images formed at 
the focal plane are feasible. Furthermore, compared to conven- 
tional pulse-echo ultrasound images, VA images generated by 
the confocal transducer do not exhibit speckle noise. 

Clinical applications of VA will demand systems based on 
linear array transducers because the ultrasound beams can be 
electronically steered over the region to be imaged. It has been 
shown through computational simulations that a typical linear 
array with 32 mm-aperture operating at 3 MHz can produce VA 
images with -6dB of lateral resolution of 1.9 mm when fo- 
cused at 5 cm in depth |16|. However, the system PSF is not 
symmetric in the imaging plane as in other VA beamforming 
schemes. The resolution is poorer in the elevation when com- 
pared to the lateral direction (roughly threefold larger). Hence, 
linear array systems (1.25 and 1.5D) may produce distorted VA 
images. Recently, a VA system based on the General Electric 



Vivid 7 ultrasound scanner with the GE 7L and lOL linear array 
transducers was designed and assessed ifTTl . This system was 
set to operate at nearly 5 MHz, with the difference-frequency 
being 51.54 kHz. When the incident beams were focused at 
25 mm in depth, the spatial resolution achieved by the system 
at -6dB was 0.573 mm laterally, 1.44 mm in elevation, and 
2.20 mm axially. 

Limitations in the spatial resolution of VA systems degrade 
the yielded images, reducing their utility to medical imaging 
applications. The need for restoration of VA images in either el- 
evation (with linear arrays) or axial directions has prompted us 
to assess the efficiency of some restoration algorithms applied 
to these images. In doing so, three deconvolution frequency do- 
main filters, namely constrained least squares (CLS), geomet- 
ric mean, and Wiener filters ||20l |21| based on the theoretical 
PSF of the system are employed to restore simulated VA im- 
ages ff 22l . In addition, a phantom composed by three stretched 
0.5 mm wires was imaged by a confocal VA system. The VA 
line spread function (LSF) was estimated from the image and 
used to restore the image through the CLS filter. Restoration 
with the theoretical LSF yielded results with no visual infor- 
mation of the wires. The main reason for this problem lies on 
the lack of the phase information in the theoretical LSF. This 
is illustrated by restoring the actual VA image with a combined 
LSF, in which the magnitude and the phase come from the the- 
oretical and the experimental LSF, respectively. In this case, 
the CLS filter restoration is capable to restore most of the infor- 
mation such as position, shape and contrast of the wires in the 
image. 

This article is organized as follows. A model for vibro- 
acoustic image formation based on a three-dimensional PSF 
composed by the product of the pressure amplitudes of the in- 
cident beams is presented in Sec. |2] The general aspects of 
the restoration filters (Wiener, CLS, and geometric mean) are 
briefly described in Sec. [3] The computational and experimen- 
tal results are explained in Sec. |4] In Sec. |5j the main conclu- 
sions of this study are outlined. 



2. Image formation 

The VA technique employs two ultrasound beams at fre- 
quencies oji and <jj2 to form an image of the region-of-interest. 
Hence the primary incident pressure upon this region is given 
by 

/7(r,0 = /^l(r)^"^"^^ + /^2(r)^"'"^^ 

where r is the position vector indicating the observation point, 
/ is the imaginary unit, t is the time, and pi and p2 are the 
complex pressure amplitudes. 

Based on the radiation force concept, a 2D image formation 
model was proposed in which the system PSF is the product of 
the incident pressure fields in the transducer focal plane fTT]. 
An extension of this model to account for depth-of-field effects 
in the VA image formation was presented in Ref. |[23ll . Accord- 
ing to these image formation models, the three-dimensional 
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PSF of a VA system is given by 

h{v) = P\{V)P2{V\ (1) 

where the symbol * means complex conjugation. It is wor- 
thy noting that the PSF of a VA system as given in Eq. ([T]) is 
a complex function. Despite early success on explaining VA 
as a radiation force imaging method, a recent study [11 J has 
shown that the difference-frequency signal used to form VA 
images depends on the nonlinear interaction of primary inci- 
dent and scattered waves at the frequencies coi and 0)2- Both 
radiation force phenomenon and the nonlinear interaction of 
the waves scale quadratically with the primary pressure ampli- 
tudes. Hence, the description of the difference-frequency pres- 
sure magnitude based upon the radiation force and the nonlinear 
interaction models are alike, but the phase information due to 
each phenomenon is different ITTII . Therefore, as the deconvo- 
lution of VA images needs to be done with the PSF (magnitude 
and phase). Thus, we can expect that the deconvolution using 
the PSF given in Eq. ([T]) may not be able to fully recover the 
visual information in VA images. 

To compute the PSF given in Eq. ([T]) one has to calculate the 
incident pressure field amplitudes pi and p2 in the linear ap- 
proximation for the VA transducer. This has been done using 
the Field II simulation program |24|. Even though the fields 
from Field II are linear the combination through equation 2 
makes the problem inherently nonlinear. The nonlinearity of 
the VA beamforming arises because of the interaction of the 
two beams, not the nonlinear propagation of the pressure from 
the transducer. The simulations were performed considering the 
transducer and operation frequencies used in the experiments 
to be described later. Hence, a two-element spherical confo- 
cal transducer with a 44 mm diameter and focus at 70 mm was 
considered. The transducer inner radius is 14.8 mm, while the 
outer ring ranges from 15.2 to 22 mm. The transducer elements 
were driven independently by sinusoidal signals at 3.075 MHz 
and 3.125 MHz. The magnitudes of the simulated PSF laterally 
and in depth are presented in Fig. [2] 

The VA image of an object can be obtained by convolving 
the function which represents the object with the system PSF. 
Accordingly, the image function is given by 1231 

g(r) = /(r)*/i(r) + ^(r), (2) 

where / is the function that represents the object, n repre- 
sents the additive Gaussian noise, and the symbol * denotes 
the spatial convolution. The additive Gaussian noise is the only 
source of noise present in VA images due to thermal fluctu- 
ations present in the electronic devices of the VA system. The 
noise level of VA images has been assessed. Thus, a white noise 
with 20 dB signal-to-noise ratio (SNR) was chosen to be added 
to all simulated images here. This level was considered by vi- 
sual inspection of the VA simulated and actual images. 

In the current arrangement of VA, real-time images cannot 
be done because the system requires four or five difference- 
frequency cycles to form a pixel. At 50 kHz, this corresponds 
to 80yus per pixel. To form images with 100 x 100 pixels and 




Figure 2: Magnitude of the PSF in the focal plane and in depth for a two- 
element confocal transducer with inner and outer radii being respectively 14.8 
and 22 mm. The transducer focal distance is 70 mm, which corresponds to zero 
in the depth scale. 



show them in 24 frames per second, one needs 4.1 yus for each 
pixel. 



3. Restoration algorithms 



Among different filtering schemes to perform image decon- 
volution, we choose the Wiener, the constrained least-squares 
(CLSs), and the geometric mean (GM) filters |20|. This choice 
is based on the following digital image processing aspects. 
Firstly, the chosen filters operate in the frequency domain, 
which is computationally less expensive than to perform the fil- 
tering in the spatial domain. Secondly, inasmuch as the VA 
image formation is based on a spatial invariant model with ad- 
ditive noise as given in Eq. ([2]), one can use inverse filter meth- 
ods for deconvolution. In this scheme, the Wiener filter re- 
moves the image degradation without division by small values, 
which can lead to numerical overflow. However, this filter re- 
duces the noise level by smoothing the image without sharp- 
ening the edges. On the other hand, the CLS filter reduces 
the image degradation and it allows the control of sharpening 
the edges through a regularization parameter. The CLS and the 
Wiener filters were used for the deconvolution of ultrasound im- 
ages |25|. The GM filter combines the features of the Wiener 
and the CLS filters depending on the adjustment of two param- 
eters. The criteria used here to set the filters parameters will be 
discussed later. 

To describe the Wiener, the CLS, and the GM filters, the im- 
age formation model in Eq. ^ is written in the frequency do- 
main as 

G(k) = F(k)/f (k) -h A^(k), (3) 

where k is the reciprocal of r, and G, //, and N are the 
Fourier transforms of the functions f,g,h, and n, respectively. 
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3.1. Wiener filter 

The Wiener filter is designed to minimize the mean squared 
error between the filtered and the ideal images. In the frequency 
domain, this filter is given by l>20J 



(4) 



where Sn and S / are the power spectrum density of the noise 
and of the ideal image, respectively. 

3.2. Constrained least-squares filter 

In the frequency domain, the constrained least-squares 
(CLSs) filter is expressed as 1201 



F = 



(5) 



where y is the regularization parameter and P is the 
Fourier Transform of the discrete Laplacian operator in three- 
dimensions. 

3.3. Geometric mean filter 

The geometric mean (GM) filter is given in the frequency 
domain by 1201 



if* 



m^ySJSf 



(6) 



where a e [0, 1] and y > 0. 



3.4. Filter parameters estimation 

In performing the VA image deconvolutions, a trial and error 
approach is adopted. In doing so, the best visual and quanti- 
tative results for the restored images are sought. The quantita- 
tive analysis is performed using three indices of digital image 
processing, namely |[26|[27l improvement in the signal-to-noise 
ratio (ISNR), mean squared error (MSB), and universal image 
quality index (UIQI). For the last two indices, the imaginary 
part of the restored image is ignored, because it approaches to 
zero after the filtering process. 

4. Results and discussion 

4. 1 . Computational simulations 

To perform the computational simulation of VA images, the 
PSF is represented by a complex matrix whose dimensions are 
256 X 256 X 512. This corresponds to a region of 32 mm lat- 
erally and 64 mm in depth. Two computational phantoms are 
devised with three and five spherical inclusions as shown in 
Fig.|3] The phantoms have dynamic range of 48 dB, i.e., a pixel 
is represented with 8 bits depth. Later, we present two com- 
puter phantoms with inclusions having up to four contrast lev- 
els. The purpose of using such phantoms is to demonstrate that 
the deconvolution of VA images can recover most information 
regarding to the inclusions morphology. Thus, the adopted dy- 
namic range, which represents inclusions with 256 : 1 contrast 



ratio, is sufficient for this proposal. In turn, the inclusions can 
be regarded as soft tissue with different mechanical properties 
compared to the background. 
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Figure 3: Phantoms used in the computational restoration, (a) Phantom 1 seen 
in the 6 mm-transverse plane, (b) Phantom 1 seen in-depth with three inclu- 
sions with 4 mm radius each and 36 dB (left), OdB (center), and 42 dB (right) 
contrasts, (c) Phantom 2 seen in the 8 mm-transverse plane, (d) Phantom 1 seen 
in-depth with two additional inclusions of 2 mm radius and 46 dB contrast. The 
dotted lines in (b) and (d) indicates the position of the images (a) and (b). 

Phantom 1 is composed by three inclusions having same ra- 
dius of 4 mm. The inclusions have contrasts of 36 dB (left), 
OdB (center), and 42 dB (right). Phantom 2 is similar to phan- 
tom 1 with two additional spheres are included with 2 mm ra- 
dius each and 46 dB contrast. The dimensions of the inclusions 
are typically found in VA images of mass legions and calcifica- 
tions. 

The VA computational images are generated by convolving 
the system PSF with the function which represents the phantom, 
as given in Eq. ([2]). A Gaussian white noise of 20 dB signal- 
to-noise ratio (SNR) was added to the generated images. The 
simulated images are presented in Fig. |4] It is not possible to 
recognize any inclusion from the axial image information. 

The simulated images are restored with the CLS, geometric 
mean and Wiener filters. The obtained results for both phan- 
toms are shown in Figs. |5] and [6] It is clear that all details in the 
axial images are recovered by the filters analyzed here. 

In Fig. [T] the profile of phantom 2 and the corresponding 
restored images at 8 mm in depth are displayed in normalized 



4 




Figure 5: Image restoration results for phantom 2 on both the 6 mm-transverse plane and in-depth. The filter parameters are y = 0.2 (CLS), and a = 0.5 and 7=1 
(geometric mean). The CLS filter yielded the best boarder preservation, while the Wiener filter yielded the most blurred image. The dotted line indicates where the 
respective transverse image is placed in depth. 



scale. The VA image profile in Fig. 7(b) shows a noisy mixture 
of two gray levels present in phantom 2. After the restoration 
process, the gray levels resemble that of phantom 2. 

The quantitative analysis is performed using the measures 



defined in Sec. 3.4 Deconvolution restoration results for phan- 
toms 1 and 2 are presented in Table [T] The values indicate good 
results for the three restoration algorithms because the ISNR 
gives 60 dB improved with respect to the raw image. Moreover, 
the MSE is lower than 0.02 and the UIQI is higher than 0.77 
(when this index values one, it indicates "perfect" restoration). 
The best restoration result was achieved by the Wiener filter 
for both phantoms. We observe that the resulting images are 
smoother than those obtained with the two other filters. In other 
words, the restored Wiener image is more blurred and the ob- 
ject borders are mixed with the background. Furthermore, the 
constrained least- squares filtered images display the best bor- 
der preservation because it is based on the Laplacian operator, 
which enhances the object borders. 

4.2. Experimental images 

To investigate how actual VA images can be restored, an ex- 
periment using a VA system is performed and the results are 
presented here. The VA system consists of the two-element 
confocal transducer described in Sec. [2] which scans a wire 
phantom (see Fig. [8]). The driving frequencies of the trans- 
ducer elements are 3.075 and 3.125 MHz and the correspond- 
ing diff'erence-frequency is 50 kHz. In turn, the wire phantom is 



composed of three stretched 0.5 mm-diameter wires distributed 
along the depth-of-field. The wires are diagonally arranged 
and spaced by 1.4 cm. The scanned area is 30 x 70 mm and 
each pixel was taken in steps of 0.1mm. This area encloses 
all depth-of-field information of the acquired images. A hy- 
drophone (ITC-6050C, ITC, Santa Barbara, CA, USA) with 
sensitivity -157 ± 3 dB re IV/yuPa is used to detect the gener- 
ated diff'erence-frequency pressure. The hydrophone is placed 
at approximately r = 20 cm way from the farthest wire to the 
transducer. This corresponds to = 41, which can be re- 
garded as the hydrophone farfield region. The active element 
of the hydrophone is a cylinder located in the middle of its con- 
ical tip. This cylinder has a height of 50 mm and diameter of 
23 mm. Its axis is aligned with the the hydrophone long axis. 
The hydrophone cannot be considered as a point detector for 
diff'erence-frequency wavelength of few centimeters (3.0 cm at 
50 kHz in water). The finite-size eff'ect of the active element (in 
a linear model) yields an output voltage proportional to the spa- 
tial average of the detected pressure over the element area plus 
a signal-independent noise | 28 1. Consequently, the hydrophone 
may add an extra information to the phase of the acquired sig- 
nal. Fig. [9] shows the experimental VA setup used. 

Two schemes are devised to restore the actual VA images. 
Firstly, we estimate the line-spread function (LSF) directly 
from the experimental image (central wire) and use it in the 
deconvolution process. In order to make the experimental LSF 
more similar to the theoretical one, we applied a spectral in- 
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Figure 6: Image restoration results for phantom 2 on both the 8 mm-transverse plane and in-depth. The filter parameters are y = 0.2 (CLS), and a = 0.5 and 7=1 
(geometric mean). The image restored with the Wiener filter is more blurred than the other results. The CLS filter yielded the best edge preservation among the 
other filters. 



Table 1: Quantitative results based on the ISNR (improvement in signal-to-noise ratio), MSB (mean squared error), and UIQI (universal image quality index) for 
phantoms 1 and 2. 



Phantom 1 



Phantom 2 



Filters 


ISNR 


MSB 


UIQI 


ISNR 


MSB 


UIQI 


Constrained least-squares 


60.13 


0.0192 


0.78 


59.91 


0.0188 


0.79 


Geometric mean 


60.74 


0.0166 


0.80 


59.83 


0.0191 


0.77 


Wiener 


61.83 


0.0129 


0.91 


61.03 


0.0145 


0.88 



version to the image, i.e., first a low-pass filtering is performed 
and then the filtered image is subtracted from the original im- 
age 1291 . Secondly, because the magnitude of both the theo- 
retical and the estimated LSF are similar, we compose the LSF 
by combining the magnitude of Eq. ([T]) and the phase infor- 
mation directly from the experimental data. Such arrangement 
was chosen because we found discrepancies between the phase 
of the estimated and the theoretical LSF based on the radia- 
tion force model. This model seems not to provide the entire 
phase information present in VA images. A more complete im- 
age formation description needs to take into account parametric 
and nonlinear scattering phenomena. 

The magnitude and phase of the estimated and the theoreti- 



cal LSFs are exhibited in Figs. lQ(a) and lQ(b) The theoretical 
LSF was obtained by the convolution of the 3D VA PSF with 



the 3D wire phantom. The phase of the experimental LSF is es- 
timated as follows. Each pixel is formed from an acquired har- 
monic sequence. This sequence discretized from the detected 
signal by the hydrophone at the diff'erence-frequency Aoj. Sup- 
pose that the time sequence which represents the /th-pixel of the 
VA image is given by 



Xi[n] = Ai cos(Aojn/fs + (pd, n = 0,1, 



(7) 



where Ai > and 0/ are the is the phase of the signal, fs = 
1 MHz is the sampling frequency, and N = 512 the number of 
sampling points. By taking the Hilbert transform 1301 of Xi[n], 
one obtains Xi[n] = A/ siniAcon / fs-\-(pi). Thus, the representation 
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(a) Phantom 2 



(b) VA image 



r 



(c) Wiener filter 



(d) CLS filter 



(e) GM filter 



Figure 7: Gray levels obtained for a transverse line at 8 mm in depth of phantom 2 . |(a)| Transverse line for phantom 2. |(b)| Transverse line for the simulated VA 
image of phantom 2. Transverse lines for the restored images using |(c) [Wiener filter,|(d)|CLS filter and|(e)|GM filter. 



of the acquired signal in complex notation is 



(8) 



Therefore, the phase of a sequence which represents the 7th- 
pixel with respect to that of the /th-pixel can be computed from 

(pj = (pi + TTIJ y\ l^^iM + ixi[n])(xj[n] - ixj[n])] . (9) 

All the pixel phases can be obtained with respect to a reference 
pixel. The amplitudes A/ and Aj can be calculated by computing 
the root mean square of the respective acquired signal. 



It can be noticed from Figs. 10(a) and 10(b) that the theoret- 
ical and the experimental LSFs have differences in their phase 
contents. The lack of information in the LSF phase leads to un- 
satisfactory restoration results with several artifacts. The miss- 
ing information in the phase is mostly related to the fact that 
the image formation model based on Eq. ([T]) does not consider 
the nonlinear interaction of the primary beams. Moreover, the 
finite- size eff'ect of the hydrophone active element was not con- 
sidered in the theoretical model. 

The result obtained by estimating the LSF from the exper- 
imental image, is presented in Fig. ll(b)| whil e the result us- 
ing the composed LSF is shown in Fig. ll(c)[ Both methods 



yielded similar restored images. Despite presenting two blurred 
wires (top left and bottom right), the restored images show as- 
pects of the imaged wires, including location, geometric shape 
and contrast. Some aspects may have contributed for the oc- 
currence of artifacts in the restored images. The deconvolution 
filters used here usually produce undesirable artifacts during the 
process of deconvolution because of the ill-posed nature of the 
image restoration problem, i.e., it may have more than one so- 
lution and the solutions depend discontinuously on the initial 
data. Therefore, perturbations in the VA image due to the noise, 
may cause undesired variations in the restoration results. This 
can be noticed by comparing the result of a complete simu- 



lated deconvolution process presented in Fig. |ll(d)[ where no 
artifacts are produced, because the whole process is controlled, 
from the image formation up until its restoration. That is, the 
same LSF is used for both VA image formation and restora- 
tion. Furthermore, it is important to note that an extra infor- 
mation may appear in the phase of the detected signal caused 
by the acoustical path diff'erence from the wire positions to the 
hydrophone. 

5. Summary and conclusions 

In this paper, we have presented the results of three restora- 
tion techniques (Wiener, CLS and geometric mean filters) ap- 
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(a) Side view. 



(b) Top view. 




16 16 

lateral (mm) 



16 16 

lateral (mm) 



Figure 4: VA images obtained by computational simulations with white noise 
of 20 dB SNR. (a) Image of phantom 1 seen in the 6 mm-trans verse plane, (b) 
Phantom 1 seen along depth, (c) Image of phantom 2 at the 8 mm-trans verse 
plane, (d) Phantom 2 shown along depth. The dotted lines in (b) and (d) indi- 
cates the position of the images (a) and (b). 




(c) View from transducer. 

Figure 9: VA experimental setup composed of a two-element confocal trans- 
ducer, a phantom with three stretched thin wires, and a hydrophone. All com- 
ponents are submerged in a water tank. 





(a) 



(b) 



Figure 10: The magnitudes of (a) the experimental and (b) the theoretical LSFs. 
The corresponding phases in radians of the center line along the depth-of-field 
are also displayed. 



Figure 8: Image of the computational representation of the wire phantom. 

plied in the restoration of VA images. The restoration algo- 
rithms based on the theoretical PSF were applied to computa- 
tionally simulated images, in which the visual information of 
the objects was lost. The performance of the algorithms were 
assessed both quantitatively and visually. The deconvolution of 
simulated images was capable of restoring the size, shape, po- 
sition and contrast of the objects present in the images for all 
filters used. 



The CLS filter was employed in restoring an actual VA image 
of three stretched wires. In this case, two LSFs were used, one 
was estimated from the image, while the other combines the 
phase of the estimated LSF with the magnitude of the theoreti- 
cal one. The wire positions and partial information about their 
shape and contrast can be recognized in the images restored 
with both experimental and combined LSFs. 

The restoration of the experimental VA image with the the- 
oretical LSF based on Eq. ([T]) yielded a poor visual restoration 
retrieve and thus was not presented here. Yet the magnitudes 
of both theoretical and experimental LSFs matched, the phases 
showed distinctly diff'erent contents. This hints that the adopted 
theoretical model for the PSF lacks the appropriate phase infor- 
mation of VA images. Additional information may come from 
the diff'erence-frequency generation based on the nonlinear in- 
teraction of the primary waves |[T0l[TTl[T2l and possibly also on 
the finite- size eff'ect of the hydrophone active element. 
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This study represents a preliminary step in understanding 
how appHcation of deconvolution filters can be used to restore 
the resolution and contrast for VA images. The optimal restora- 
tion filter may depend on the object being imaged with to im- 
prove resolution and contrast. Future studies would be neces- 
sary to understand and validate how diff'erent filters and mod- 
eling of the PSF would aff'ect the restoration of VA images of 
simple phantoms and soft tissues. 
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Figure 11: Image restoration using the CLS filter with y = 0.01. The raw VA 
image is shown in (a). The restoration results are presented for the experimental 
(b) and the composed (c) LSFs. The deconvolution for the computer generated 
VA image is shown in (d). 
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